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DFT(B3LYP, B3PW91) level calculations in conjunction with three different of basis sets have 
been used to investigate the variations in the bond lengths, dipole moment and rotational constants, 
IR frequencies, IR intensities and rotational invariants of CICCCN. The nuclear quadrupole con- 
stants of chlorine and nitrogen of CICCCN have been calculated on the experimental r^ structure 
as well as on the B3PW91/6-311-I— |-g(d,p) optimized geometry and are found to be within the scale 
length of the experimental uncertainty. The slopes and intercepts obtained from the regression 
analysis between the B3LYP/6-311-I— |-g(d,p) level calculated and experimental Bo values of 
CICCCN were used to calculate the reasonable values of rotational constants of all the rare isotopic 
species of CICCCN having standard deviation ±0.048 MHz. All the spectroscopic parameters 
obtained from DFT calculation shows satisfactory agreement with the available experimental data. 



I. INTRODUCTION 



Quantum chemical calculation and spectroscopic char- 
acterization of organic compounds, free radicals, radical 
anions etc. have found considerable amount of interest 
in recent yeaiai'^'^. Accurate and efficient calculation 
of spectroscopic constants for a wide range of molecu- 
lar systems employs readily available methods and ba- 
sis sets^. The DFT method has been demonstrated to 
have great accuracy in reproducing the experimental val- 
ues of quadrupole hyperfine coupling constants^, molecu- 
lar structural properties^, rotational constantsi*^, IR fre- 
quencies, IR intensities and rotational invariantsSiiSiii 
etc. to within 8-10%. DFT studies have been car- 
ried out for many compounds like BrCCCNiSi, FCgHi^, 
ClCNiii etc. and benzene derivatives like o-benzyne^, 
S-chlorobenzonitrilsi^ etc. to calculate their accurate 
spectroscopic parameters. Highest number of the poly- 
atomic molecules identified in space are mostly carbon 
chainsiSiii, the longest chain identified so far is HCn?^. 
Identification and characterization of these species could 
have been made by their previously reported spectro- 
scopic parameters. Spectroscopic studies of the related 
series of molecules containing halogens have been exten- 
sively studiediSiSSiSiiS^. Recently, we have investigated 
the variations of the spectroscopic parameters of bromo- 
cyanoacetylene (BrCCCN) at the HF-SCF and DFT lev- 
els in conjunction with a variety of basis sets^^. It was 
found from the the investigation that DFT-B3LYP/6- 
311-1— |-g(d,p) level calculation is satisfactory in the pre- 
diction of bond lengths and rotational constants etc. of 
BrCCCN when compared against their experimental val- 
ues. However, to my knowledge, a detailed quantum 
chemical calculations and molecular properties of chloro- 
cyanoacetylene (CICCCN) is very limited. Analysis of 
the microwave spectra of CICCCN has been reported by 
T. BjorvattenS^. Infrared studies of CICCCN have been 
reported by S.J. Cyvin et.al^i. Later on, P. Klaboe et. 



4000 cm^i. 

In this study, the results of DFT (B3LYP, B3PW91) 
level calculations of optimized molecular geometry, rota- 
tional constants, quadrupole coupling constants, IR fun- 
damental frequencies, IR intensities and rotational in- 
variants of CICCCN are compared against their available 
experimental data and reported. Linear regression anal- 
ysis between the B3LYP/6-311-I— l-g(d,p) level calculated 
and experimental Bq values have been made for a reason- 
able prediction of rotational constants of all the other 22 
rare isotopomers of CICCCN. Due to the lack of quan- 
titative data on the IR intensities in the literature, the 
rotational invariants^LSSiSS of CICCCN have been pre- 
dicted and discussed here may be helpful in future for the 
experimental IR spectroscopist in interpreting IR inten- 
sities of this molecule. Satisfactory agreements between 
the calculated and available experimental values of spec- 
troscopic parameters have been found at the B3LYP level 
and compared with the corresponding values calculated 
at B3PW91, RHF and MP2 levels of theory In addition, 
a test of the optimization followed by a frequency calcula- 
tion at the MP2 level in conjunction with a medium size 
basis set 6-31G or/and 6-311++g(d,p) predicts two nega- 
tive frequencies (saddle point of oder two) corresponding 
to the -C=C-C1 doubly degenerate bending mode (1^7) 
and thus it was difficult to calculate the harmonic fre- 
quencies of these modes under Cu point group symmetry. 

II. COMPUTATIONAL METHODS AND 
CALCULATIONS 

Geometry optimization and quantum chemical calcula- 
tions were carried out at the restricted HF-SCF method, 
MP2 method and Density Functional Theoretical (DFT) 
method under C^ point group symmetry. The hybrid 
HF/DFT methods used were Becke's three-parameter 
method^, with Lee- Yang-Parr correlation (B3LYP), 
Becke's three parameter exchange with Perdew-Wang 
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Five different basis sets used during the analysis were 
Dunnings correlation consistent polarized valence double 
and triple zeta basis sets aug-cc-PVNZ (N = D, T)^** 
augmented with a d function, Ahlriches TZV(3df,2p)'^^ 
and triple split zeta qualities 6-311H — |-g(d,p) and 6- 
311+g(df,pd) augmented with p, d and f functions. 
Tight convergence in stead of default convergence 
criteria was used with ultrafine integration grid for the 
calculation of normal mode frequencies. All calculations 
were carried out by the GaussianOSW suite program 
package2&. 

The physical quantities pa, (3a and Xa sjie invariant 
with respect to the rotation of space-fixed coordinate 
axes which are known as generalized atomic polar tensor 
charge (GAPT)ii, atomic anisotropy and Kings atomic 
effective charge of ath atom respectively and are useful 
for interpreting infrared intensities. The components of 
the (3x3) atomic polar tensor matrix pjf'' are defined 
as the first derivatives of the components of molecular 
dipole moment with respect to the atomic Cartesian 
displacement co-ordinates of each atom a (=C1, C and 
N) of CICCCN given by: 



where L~^,\J and B are matrices^i used for molecular 
vibrational analysis and Pq contains the dipole moment 
derivatives with respect to the normal co-ordinates, 
which are proportional to the experimental infrared 
intensities^^. Pp is the rotation plus translation polar 
tensor whose elements for a neutral molecule propor- 
tional to and p = BX, where n and I are dipole 
moment and moment of inertia of of the molecule 
respectively. Atomic polar tensor components were 
directly calculated by Gaussian03 package. 

On the other hand, nuclear quadrupole coupling 
constants (NQCCs) of chlorine and nitrogen of CICCCN 
are calculated by using the electric field gradients 
(EFGs). The elements of the NQCC tensors Xij 
related to those of the EFG tensors qij by : 

Xij = {eQ/h)q,j (8) 

where e is the fundamental electronic charge, h is 
Plank's constant and i,j — a,b,c are the principal axes 
of the inertia tensor. 
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The generalized atomic polar tensor charges (GAPTs) 
are nothing but mean dipole derivatives^^iS^ given by: 
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for the neutral molecule CICCCN. The quantity 
known as the effective charge of the ath atom, defined 
by King and co-workers^SiSS as: 



xi-l[Tr.iP'pPT^)]-pl + m 



(4) 



where 
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The molecular polar tensor of CICCCN is a juxtapo- 
sition of the (3x3) atomic polar tensors given by: 



Px ^ {p^^\p^\p^\p^\pn 

and is calculated from 

Px = PqL-'UB 



(6) 



III. RESULTS AND DISCUSSION 

Optimization of geometry for each molecule at the 
restricted HF-SCF and DFT levels were performed in 
the ground state by single point energy calculation. 
In each case, stationary points were found and linear 
geometry of CICCCN was confirmed. Finally, a fre- 
quency calculation following each optimization have 
been performed in order to check for the existence of a 
true minimum and to confirm an equilibrium structure. 
A test of geometry optimization followed by a frequency 
calculation under point group symmetry at the MP2 
level in conjunction with 6-31g or/and 6-311++g(d,p) 
basis sets succeeded in locating a stationary point but 
with two imaginary frequencies (IMAG = 2) correspond- 
ing to the -C=C-C1 doubly degenerate bending mode. 
This confirms that the stationary point located by the 
MP2 level of theory is not a true minimum rather a 
saddle point of order two on the potential energy surface 
of CICCCN which may cause a large deformation to 
the geometry of the molecule. Similar features have 
also been reported elsewhereSS. So the scan of the 
potential energy surface should be emphasized in order 
to characterize the the nature of the saddle point on 
the potential energy surface of CICCCN and to check 
any possible existance of conical intersection between 
the ground and first excited state. However, the DFT 
and HF-SCF level calculations were not able predict any 
imaginary frequencies implying that the stationary point 
is located at the global minimum of the potential energy 
hyper-surface. The results of our calculation for bond 
lengths of CICCCN are summarized in Table 1. The 
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TABLE L Comparison of the molecular optimized geometry (in A), dipole moment, rotational constant and total energy (a.u.) 
of CICCCN (calculated by various methods in conjunction with basis sets of increasing size) with the experimental average Ts 
values^'' . 



Method 


basis 


Cl-C 




C-C 




C1..N 




B/MHz 




AO. a 


Energy 


HP 


6-311++g(d,p) 


1.6359 


1.1809 


1.3855 


1.1355 


5.3378 


3.82 


1394.78 


-81.89 


-3.81 


-627.5026 




aug-cc-pVDZ 


1.6421 


1.1881 


1.3914 


1.1372 


5.3588 


3.85 


1381.38 


-80.34 


-3.61 


-627.4784 




aug-cc-PVTZ 


1.6325 


1.1780 


1.3852 


1.1277 


5.3234 


3.86 


1399.53 


-82.34 


-4.30 


-627.5349 


MPa 


6-311++g(d,p) 


1.6294 


1.2238 


1.3705 


1.1826 


5.4063 


4.12 


1362.79 






-628.1965 


B3LYP 


6-311++g(d,p) 


1.6322 


1.2064 


1.3647 


1.1592 


5.3625 


4.17 


1382.60 


-72.82 


-3.48 


-629.2369 




aug-cc-pVDZ 


1.6402 


1.2159 


1.3721 


1.1670 


5.3952 


4.17 


1366.05 


-77.31 


-2.10 


-629.2044 




aug-cc-pVTZ 


1.6276 


1.2039 


1.3643 


1.1562 


5.3520 


4.14 


1380.00 


-76.10 


-3.64 


-629.2646 


B3PW91 


6-311++g(d,p) 


1.6232 


1.2073 


1.3632 


1.1595 


5.3532 


4.19 


1387.85 


-74.93 


-3.38 


-629.1045 




aug-cc-pVDZ 


1.6312 


1.2160 


1.3700 


1.1671 


5.3843 


4.13 


1371.95 


-76.05 


-2.84 


-629.0762 




aug-cc-pVTZ 


1.6196 


1.2049 


1.3630 


1.1567 


5.3442 


4.16 


1392.34 


-74.16 


-3.54 


-629.1332 


Expt. 




1.6245 


1.2090 


1.3690 


1.1602 


5.3627 


3.38 


1382.328(2) 









TABLE II: Comparison of calculated" and experimental— values of quadrupole coupling constants Xaa{= Xzz) MHz) of 
^^Cl, ■^^Cl and '^''N of CICCCN. The values given in the parentheses are directly calculated without using the scaling factorai. 



Atoms 




XB3LYP 


XB3PW91 


XExpt. 


35 Q 


-79.57(-80.39) 


-80.67(-81.48) 


-79.53(-80.35) 


-75±4 


3^C1 


-62.94(-63.50) 


-63.58(64.14) 


-62.68(-63.25) 


-62±3 


14n(35c1) 


-4.29 


-4.30 


-4.31 




i''N(^^Cl) 


-4.30 









" See text for discussion. 

^ Refi-2<. Actual bond distances based on the rs co-ordinates are reported to be 1.6233(31) (1.6256(32)), 1.2086(42) 
(1.2093(43)), 1.3700(23) (1.3680(21)), 1.1607(10) (1.1596(10)) for Cl-C, C=C, C-C, C=N of 35 (and 37) species of CICCCN 

respectively. 



level calculations highly overestimates Cl-C and C-C 
bond lengths and underestimates C=C and C=N bond 
lengths, whereas these values calculated at the MP 2 level 
overestimates all these bond lengths w^hen compared 
against their experimental values. However, the results 
of DFT-B3LYP, B3PW91 method in conjunction with 
6-311-1— |-g(d,p) basis set are in rather pleasing agreement 
with the experimental Ts values. Here, We note a slight 
overestimation of Cl-C bond length and an underestima- 
tion of C=C and C=N bond lengths, differences between 
calculated values and the average rs values in bond 
lengths for CICCCN are 0.0077 A, -0.0026 A, -0.0043 A 
and -0.001 A for Cl-C, C=C, C-C and C=N respectively 
at the B3LYP level, whereas a slight underestimation of 
all these bond distances being noticed at the B3PW91 
level of theory, differences are -0.0013 A, -0.0017 A, 
-0.0058 A, -0.0007 A for Cl-C, C=C, C-C and C=N 
respectively. A comparison of bond lengths calculated 
at various levels, summarized in Table 1, indicate 
that convergence has not achieved on improving the 
size of the basis sets 6-311-1— |-g(d,p) to aug-cc-PVTZ. 
However, a better comparison is made on the over-all 
(end-over-end) length , rci...N, of the molecule which is 
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FIG. 1: Comparison between the calculated (B3LYP/6- 
311++g(d,p)) and experimental overall geometry of CICCCN. 



B3LYP level in conjunction with 6-311++g(d,p) basis 
set is calculated to be 5.3625 A which differs from 
the average experimental value by a factor of 0.0002 
A . Thus, it is to say that the over-all r^-structure 
is a good approximation to the over-all bond length 
at the correlated level B3LYP/6-311-f-|-g(d,p) within, 
say, 0.001 A , which is mostly under the statistical 
uncertainties. 

Interestingly, our B3LYP/6-311++g(d,p) level 
optimized geometry leads to satisfactory values of 
rotational constants of CICCCN shown in Table 3. 
The experimental values are also given for compari- 
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TABLE III: Comparison of the rotational constants of all the 
isotopomers of CICCCN calculated at B3LYP/6-311++g(d,p) 
level against their experimental values^"^ . 



Isotopic Species 


Calc. 


Scaled" 


Expt. 


CICCCN 


1382.602 


1382.339 


1382.328(2) 


^^Cl^^CCCN 


1381.751 


1381.488 


1381.454(2) 


^^CIC^^CCN 


1380.607 


1380.346 


1380.357(2) 


^'^CICC^^CN 


1366.338 


1366.090 


1366.054(2) 


^^CICCC^^ N 


1344.232 


1344.005 


1343.913(2) 


^^Cl^^CCC^^N 


1343.295 


1343.068 




^^CIC^^CC^^N 


1342.535 


1342.309 




^^CICC^^C^^N 


1329.395 


1329.181 




^^Cl^^C^^CCN 


1379.727 


1379.466 




■^^Cr^C^CC^N 


1341.571 


1341.346 




^^Cl^^CC'^^CN 


1365.420 


1365.173 




^^Cl^'^CC^^C^^N 


1328.391 


1328.178 




^^CIC'^^C^^CN 


1364.517 


1364.271 




^^CIC^^C^^C^^N 


1327.846 


1327.634 




^^Cl^'^C^^C^'^CN 


1363.571 


1363.325 




^■^ Cl^^ C^^ C'^^ C^^ N 


1326.817 


1326.790 




^"^CICCCN 


1350.544 


1350.311 


1350.360(2) 


^''Cf^CCCN 


1349.888 


1349.655 


1349.692(2) 


^"^ CIC^^CCN 


1348.379 


1348.148 


1348.202(2) 


^'^CICC^^CN 


1334.300 


1334.082 


1334.089(2) 


^''CICCC^^N 


1312.831 


1312.633 


1312.603(2) 


^"^Cl^^CCC^'^N 


1312.099 


1311.901 




^''CIC^^CC^^'N 


1310.977 


1310.780 




^''CICC^^C^^'N 


1298.001 


1297.827 




^"^Cl^^C^^CCN 


1347.697 


1347.466 




^"'Cl^'^C^^CC^'^N 


1310.220 


1310.024 




^''Cl^^CC^^CN 


1333.585 


1333.367 




^''Cl^^CC^^C^^N 


1297.209 


1297.025 




^''CIC^^C^^CN 


1332.317 


1332.108 




^''CIC^^C^^C^^N 


1296.302 


1296.119 




^'^Cl^^C^^C^^CN 


1331.576 


1331.360 




^'^Cl^^C^^C^^C^^N 


1295.486 


1295.304 




" Scaled using the slop 


es and intercepts obtained from the 


linear 


regression 


analysis. 





%) than the experimental Bq values of all the 10 
different isotopomers, whereas these values calculated 
in conjunction with other basis sets are far off (sum- 
marized in Table 1). Thus, the agreement between the 
theoretical prediction and experimental observation of 
rotational constants clearly indicates that the B3LYP /6- 
311++g(d,p) optimized average geometry, shown in 
Fig 1, is realistic and advocates for a decent predictive 
power of the computational procedure applied to CIC- 
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1305 1320 1335 1350 1365 1380 
Calc. 

FIG. 2: Comparison between the calculated (B3LYP/6- 
311+- |-g(d,p)) and experimental rotational constants of CIC- 
CCN. The slopes and intercepts obtained from the linear re- 
gression analysis are 0.99907 and 1.02283 respectively. 



analysis, shown in Figures 2, between the calculated and 
experimental Bo values, rotational constants of other 
22 rare isotopomers of CICCCN has been predicted 
accurately having standard deviation ±0.048 MHz. A 
closer look at the Table 1 indicates that the HF-SCF 
wave functions predict dipole moments of CICCCN close 
to its experimental value, deviation being 0.4 D where 
as the DFT wave functions overestimates it by a factor 
of 0.72 D irrespective of the basis sets used. Zero point 
vibrational effects are not taken into consideration in all 
these calculations. 

The nuclear quadrupole coupling constants (NQCCs) 
of 35C1, 37C1 and of CICCCN, NQCC being propor- 
tional to the the electric field gradient (EFG), are calcu- 
lated by using equation (8). Calculation of EFGs being 
made on the rs-structures^S, as well as on the B3LYP 
and B3PW91/6-311++g(d,p) optimized geometries of 
CICCCN. The methods, basis sets and procedures^ used 
during the calculation were B3PW91/6-311-|-g(df,pd) 
for nitrogen and BlLYP/TZV(3df,2p) for chlorine 
respectively. The coefficients (eQ/h)eff. in equation (8), 
obtained from the linear regression analysis between 
the calculated EFGs and experimental NQCCs for a 
series of molecules containing chlorine and nitrogen, are 
reported^ to be 4.5586(40) MHz/a.u., -19.185 MHz/a.u. 
and -15.120 MHz/a.u. for ^^n, 35^1 and ^^Cl respec- 
tively. The nuclear quadrupole coupling constants of 
3^C1, ■^''Cl and ^''N are calculated on the structure is 
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geometries are also included in Table 2 for comparison. 
As can be seen, the values of NQCCs deviates from the 
experimental values by a factor of 4.57 MHz (6.09 %) 
and 0.94 MHz (1.5 %) for ^^C\ and ^^Cl respectively. 
These deviations are significantly due to the largest 
error ±4 MHz and ±3 MHz present in the measurement 
of the quadrupole hyperfine structures of '^^Cl and 
^"^Cl of CICCCN, though the calculated values are 
within the scale length of experimental uncertainties. 
A comparison of NQCCs of Chlorine calculated on the 
B3PW91/6-311++g(d,p) geometry as well as on the 
rs structure reveals an almost similar value, which is 
basically due to the similar C-Cl bond length, whereas 
its value calculated on the B3LYP/6-311++g(d,p) 
optimized geometry is a bit off as expected since NQCC 
varies linearly with C-Cl bond length. A comparison 
with ClCNi^ shows an underestimation in the value 
of NQCC of CI, which being supported by a decrease 
of C-Cl bond length, the decrease being within 5%. 
On the other hand, the nuclear quadrupole coupling 
constant of ^'*N of CICCCN is also summarized in Table 
2. As can be seen, the values of NQCCs of ^''N are very 
similar to each other since the r^ structure as well as 
the optimized B3LYP, B3PW91/6-311++g(d,p) C=N 
structures are predicted to be very similar. Finally, 
comparing FCCCNl, HCCCN and DCCCN^ with that 
of CICCCN reveals a similar C=N structure resulting 
a similar value of NQCC. Thus, the accurate values of 
NQCCs of 35C1, 37ci and are those derived from 
the r^ structures. Variations of NQCCs of '^'^Cl and N 
of CICCCN calculated directly by various levels are also 
summarized in Table 1. 

CICCCN has, in principle, ten fundamental modes. 
These modes under the point group symmetry are 
distributed among four stretching vibrations of E'^ 
species and six bending vibrations of H species, out of 
which three H species are doubly degenerate. Thus, 
the vibrational assignments of CICCCN reported by the 
co-authors2ii2Si2& comprises seven fundamental modes 
J^i , 1^2, J^5, 1^6 and vj which is in consistent with 

the HF and DFT level calculations. All the normal 
modes are calculated to be Raman active. Harmonic 
frequencies have been calculated at the RHF and DFT 
levels in conjunction with three different basis sets with 
increasing size in solution (benzene) as well as in vapour 
phase are summarized in Table 5. Experimental values 
are also included for comparison. All the predicted vi- 
brational spectra have no imaginary frequency, implying 
that the optimized geometry is locating at the global 
minimum of the potential energy hyper-surface for both 
the methods. Frequencies calculated at the Hartree-Fock 
level contain known systematic errors due to the lack 
of electron correlation and the choice of the basis sets 
used, resulting an overestimation of 10-12 %. Therefore, 
it is usual to scale stretching frequencies predicted at 



values to a high degree accuracy for a wide range of 
systems'*^. However, the improvement of the quality of 
the unsealed normal mode frequencies is significant in 
going from the non-correlated to the correlated level of 
theory and enlarging the size of the basis sets at the 
cost of its computational expanse. A comparison of the 
experimental values against the resulting IR normal 
mode frequencies calculated at the B3LYP level leads 
to a largest over-all error with 6-311-|-+g(d,p) basis 
set is 32 cm~^, with the aug-cc-pVDZ basis set the 
over-all error is 16 cm^^ whereas with aug-cc-pVTZ 
the over-all error is 18 cm~^. An unexpected change 
of the fundamental frequency has been noticed both 
at the B3LYP and B3PW91 levels of calculation when 
6-311++g(d,p) results being compared against the 
aug-cc-pVNZ (N — D, T) results indicating that this 
mode is sensitive to the size of the basis sets used. Thus, 
the size of the aug-cc-pVDZ basis set is somewhat more 
precise for the prediction of normal mode frequencies 
of CICCCN at the DFT levels of theory. The values of 
vii 1^3, 1^4 and computed at B3LYP/aug-cc-pVDZ 
are overestimated from their respective experimental 
values by 85.1 cm~^, 49.8 cm~^, 25.1 cm~^, 13.7 cm~^, 
13.2 cm~^ , whereas vq and i^y are underestimated by a 
factor of 21.4 cm^^ and 7.4 cm^^. These discrepancies 
can be corrected either by computing an-harmonic 
force constants or introducing a scaled field or directly 
scaling the calculated wave-numbers by taking the ratio 
between the calculated and observed frequencies for 
a particular type of motion. B3LYP scaling factor 
for stretching modes are all close to 0.965 published 
elsewhere^. Similarly, the normal mode frequencies 
calculated at the DFT-B3PW91 level of theory slightly 
overestimated from the DFT-B3LYP values, differences, 
for example, being 2.9, 2, 1.3, 8.8, 15.6, 9.3 and 13.6 
cm~^ for 1^7. ..and vi respectively in conjunction with 
aug-cc-pVDZ basis set, whereas this overestimation is 
highly overestimated for stretching and bending modes 
computed at the HF level of theory, differences being 
20.9, 86.2, 84.7, 23.6, 36.2, 235.1 and 269.7 cm"! respec- 
tively. A comparison of the normal mode frequencies 
calculated at the B3LYP/aug-cc-pVDZ in solution (ben- 
zene) with those of gas phase values indicates a shift of 
C=N stretching frequency towards the high wavelength 
region as expected whereas for other modes the shift is 
rather small which being in good agreement with the 
experimental shifts. On the other hand, a frequency 
calculation at the MP2 level in conjunction with a 
medium size basis set 6-311-|--|-g(d,p), presented in Table 
5, predicts two negative frequencies (saddle point of oder 
two) corresponding to the -C=C-C1 doubly degenerate 
bending mode (1/7) and thus it was difficult to calcu- 
late the harmonic frequencies of these modes at this level. 

The infrared band intensities li, on the other hand, of 
CICCCN corresponding to the seven vibrational modes 
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TABLE IV: RHF, B3LYP and B3PW91 level (in conjunction with 6-311++g(d,p) (A), aug-cc-PVDZ (B), aug-cc-PVTZ (C)) 
calculated atomic mean dipole derivatives (pa), tensor anisotropies {/3a), charge undeformability {pa/ Pa) and effective charges 
{Xa), Mullikan Charges {^' qa) of ChC2Ci,CAN5 in the unit of e. 

Invariants RHF B3LYP B3PW91 





A 


B 


C 


A 


B 


C 


A 


B 


C 


Pci 


-0.144 


-0.128 


-0.126 


-0.158 


-0.148 


-0.145 


-0.165 


-0.154 


-0.150 


Pc 


0.398 


0.383 


0.385 


0.398 


0.394 


0.395 


0.410 


0.4021 


0.403 


pc 


-0.259 


-0.266 


-0.276 


-0.228 


-0.243 


-0.252 


-0.231 


-0.247 


-0.257 


Pc 


0.395 


0.408 


0.418 


0.341 


0.345 


0.359 


0.333 


0.342 


0.354 


pN 


-0.390 


-0.397 


-0.401 


-0.353 


-0.348 


-0.357 


-0.347 


-0.340 


-0.350 


I3ci 


0.753 


0.725 


0.725 


0.926 


0.873 


0.864 


0.956 


0.896 


0.886 


Pc 


1.246 


1.229 


1.231 


1.430 


1.338 


1.329 


1.465 


1.361 


1.352 


Pc 


0.835 


0.813 


0.816 


0.890 


0.817 


0.809 


0.903 


0.819 


0.814 


Pc 


0.556 


0.487 


0.478 


0.558 


0.490 


0.475 


0.562 


0.483 


0.467 


Pn 


0.214 


0.179 


0.168 


0.173 


0.139 


0.131 


0.168 


0.129 


0.119 


Pci/P 


-0.191 


-0.176 


-0.174 


-0.171 


-0.169 


-0.168 


-0.173 


-0.172 


-0.169 


Pc/P 


0.319 


0.312 


0.313 


0.278 


0.294 


0.297 


0.280 


0.295 


0.298 


Pc/P 


-0.310 


-0.327 


-0.338 


-0.256 


-0.297 


-0.311 


-0.256 


-0.302 


-0.316 


Pc/P 


0.710 


0.838 


0.874 


0.611 


0.704 


0.756 


0.592 


0.708 


0.758 


Pn/P 


-1.822 


-2.218 


-2.387 


-2.040 


-2.504 


-2.725 


-2.005 


-2.636 


-2.941 


\xci\ 


0.383 


0.365 


0.364 


0.464 


0.437 


0.432 


0.480 


0.449 


0.444 


\xc\ 


0.710 


0.694 


0.697 


0.783 


0.744 


0.741 


0.804 


0.757 


0.754 


\xc\ 


0.471 


0.467 


0.473 


0.477 


0.455 


0.457 


0.483 


0.458 


0.462 


\xc\ 


0.474 


0.469 


0.474 


0.430 


0.416 


0.423 


0.425 


0.410 


0.417 


IXivl 


0.403 


0.406 


0.409 


0.362 


0.354 


0.362 


0.355 


0.348 


0.354 


qci 


0.387 


-0.208 


-0.063 


0.463 


-0.030 


0.042 


0.510 


-0.137 


-0.003 


qc 


0.082 


-0.898 


-0.901 


0.035 


-0.661 


-0.877 


0.101 


-0.678 


-0.820 


qc 


1.954 


1.838 


1.985 


1.648 


1.155 


1.656 


1.749 


1.367 


1.529 


qc 


-2.139 


-0.387 


-0.632 


-1.969 


-0.214 


-0.447 


-2.154 


-0.214 


-0.252 


qN 


-0.283 


-0.344 


-0.389 


-0.177 


-0.249 


-0.374 


-0.207 


-0.338 


-0.454 



very strong, strong, very week, weak, strong, medium 
strong, and very weak respectively. The calculated in- 
tensities are summarized in Table 7. As can be seen, the 
most intense absorption infrared band (Ii) corresponds 
to the C=N normal mode as expected irrespective of 
the methods and basis sets used, variation in intensity is 
within 30 km/mole between the two (HF-SCF and DFT) 
theoretical methods. For the band I2, the differences 
between the RHF and B3LYP values in intensity are 
within 20 km/mole. Out of the three doubly degenerate 
bending modes, intensities of /g and J5 corresponding 
to the normal modes z^s and are calculated to be of 
equal in magnitude which were experimentally assigned 
to be medium strong and strong respectively. As shown 
in Table 6, the main discrepancy in intensity is found for 
the =C-C= stretching band I3. For this band, variations 
in intensities are within 6-13 km/mole between the two 
theoretical methods. A comparison of IR intensities 
calculated (for example, at the B3LYP/aug-cc-pVDZ 



indicates a significant enhancement of C=N stretching 
intensity as expected, since intensities measured in solu- 
tions somewhat larger than the gas phase intensitiesiS, 
whereas for other bands the change is very less. Thus, 
the agreement between the calculated and the experi- 
mental gauss treatment in intensities, is far from being 
quantitative. Since there is no quantitative experimental 
data on IR intensity measurements for all these bands, 
no satisfactory interpretation could have possible. How- 
ever, a direct calculation of the dipole moment derivative 
with respect to the normal co-ordinates can serve the 
purpose. For, one makes the use of APT analysis that 
permits the interpretation of IR intensities by means 

of mean atomic charges (pa), the mass- weighted square 

2 

effective charges (^) , where m^ and Xa are the mass 
and effective charge of each atom a of CICCCN, since 
all these invariants are related to each other by the well 
known relation called as the G-sum rule^ given by: 
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FIG. 3: Normal modes (summarized in Table 5) of CICCCN. 
Atomic numbering is shown in Figure 1. The modes corre- 
sponding to ;/5 — v-j are repeated twice since these are doubly 
degenerate. The actual directions of vibrations of the doubly 
degenerate bending modes were out of (shown as up) and into 
(shown as down) the paper which being modified here for a 
clear visualization of these modes. 



it permit a direct calculation of atomic polar charges. 

The mean atomic charges as well as the other 
rotational invariants derived from the atomic polar 
tensors (APTs) using equations (1) to (5) for all the five 
atoms of the CICCCN are summarized in Table 4. The 
atomic polar charges Pa represent the redistribution of 
the electric charge density around each atom in such a 
way that their sum over all atoms equal to to zero in the 
molecular bonding environment of CICCCN at its op- 
timized equilibrium geometry. An important difference 
between APT charges and the MuUikan charges^^i^ 
is that the basis-set dependence of the former arises 
only from the fact that the basis set can be incomplete; 
hence, as the basis set approaches completeness, the 
APT charges approach a well-defined limit''^. A test 
with CICCCN, shown in Table 4, clearly indicates that 
GAPTs are almost invariant with respect to the basis 
sets. The polar tensor elements are calculated to be all 
negative for chlorine and nitrogen atoms as expected 
since these are the electronegative elements, contrast to 
the carbon atom (C4) attached to it through a triple 
bond character which is highly positive, resulting the 
mean dipole moment of nitrogen as negative and that of 
carbon as positive, each of which can be interpreted as 
carrying slightly more than -0.35 e and -1-0.35 e. This, 
in turn, results in an unequal distribution of charges 



almost equal amount of positive charge on the other 
carbon atom C2, leading to an interpretation that 
the carbon atoms C3 and C2 carries atomic charges 
slightly more than -t-0.25e and -0.25e respectively. The 
average values of atomic polar charges calculated by 
all the methods are -0.146e, 0.396e, -0.251e, 0.366e and 
-0.365e having standard deviations ±0.013e, iO.OOle, 
±0.0.016e, ±0.032e and ±0.024e for CI, C, C, C, N of 
CICCCN respectively which indicates an overestimation 
of standard deviations of C4 and N that arises mainly 
due to Y>zz component of the atomic polar tensors of 
C and N. It is noted worthy that all other invariant 
quantities are within O.le whereas a larger variation in 
the value of anisotropy parameter (3 as well as in the 
value of undeformability of charge [P / (3) on the atom a 
of CICCCN (as shown in Table 4) has been noticed. 

The Mulliken atomic charges, on the other hand, of 
chlorine and carbon attached each other by means of 
a single bond are predicted to be positive in conjunc- 
tion with 6-311-1— |-g(d,p) basis set at all the methods, 
whereas these are predicted to be negative in conjunction 
with aug-cc-pVNZ (N — D, T) basis sets. Similarly, 
the charge of carbon atom attached to the nitrogen 
atom through a triple bond character is predicted to be 
negative irrespective of the methods and basis sets used. 
However, it successfully predicted a real sign to the 
nitrogen atom as expected. As can be seen, variations 
of these charges with respect to all the methods as a 
function of basis set choice is rather large, variations 
being -0.208e to 0.51e, -0.898e to O.lOle, 1.155e to 
1.985e, -0.387e to -2.154e and -0.117e to 0.454e for CI, 
C, C, C and N of CICCCN respectively, whereas APT 
charges are less variable ireespective of the methods and 
basis sets used. A comparison of APT charges with the 
Mulliken atomic charges indicates that the APT charges 
are highly reliable than the Mulliken atomic charges 
and the Mulliken population analysis of atomic charge 
densities of CICCCN does not refiect the physical and 
chemical characteristics of this system under consid- 
eration. Finally, a comparison of the Mulliken atomic 
charges with the effective charges (%), shown in Table 4, 
reveals an increase of Chlorine(Cl), C2, C4 and A^s APT 
charges as diminishes, which leads to an interpre- 
tation that X is dominated by the stretching intensities^. 



IV. CONCLUSION 

All the spectroscopic constants were calculated at the 
restricted HF-SCF as weU as DFT (B3LYP, B3PW91) 
levels in conjunction with a variety of basis sets. Satisfac- 
tory agreements between the B3LYP/6-311-|--|-g(d,p) and 
experimental values of rotational constants were found 
for CICCCN. Over-all, rci...N, bond distance is a good 
approximation to the over-all bond distance at the corre- 
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TABLE V: Comparison of the normal mode frequencies" (in cm~^) of CICCCN calculated at all the three different levels 
of theory with their corresponding experimental values. Ifere A, B, and C stands for 6-311++g(d,p), aug-cc-PVDZ and 
aug-cc-PVTZ basis sets respectively. 



fi Species 


Motion'' 




RHF 






B3LYP 






B3PW91 




MPa 


Expt.^Expt.'' 








A 


B 


C 


A 


B 


C 


A 


B 


C 


A 






1 


E+ 


~C =N St. 


2652.9 


265.9 


2644.3 


2384.4 


2382.2 


2375.9 


2397.3 


2395.8 


2386.7 


2263.7 


2297 


2283 


2 


E+ 


-C=C St. 


2477.0 


2478.9 


2474.6 


2248.3 


2243.8 


2249.2 


2256.8 


2253.1 


2256.0 


2087.6 


2194 


2195 


3 


E+ 


eeC-C= St. 


1152.6 


1154.4 


1148.4 


1116.5 


1118.2 


1115.1 


1132.6 


1133.8 


1129.5 


1112.6 


1093 


1099 


4 


E+ 


=C-C1 St. 


563.4 


564.4 


563.0 


540.1 


540.8 


541.3 


549.1 


549.6 


549.3 


538.2 


527 


528 


5 


n 


eeC-C=N bd. 


582.8 


580.9 


615.6 


507.3 


496.2 


535.2 


511.9 


497.5 


536.8 


412.6 


483 


482 


6 


n 


=C-C=C- bd. 


379.1 


397.8 


423.5 


268.2 


311.6 


361.6 


268.0 


313.6 


364.8 


129.9 


333 


332 


7 


n 


-C=C-C1 bd. 


150.0 


142.4 


150.5 


135.7 


121.5 


137.0 


135.4 


123.4 


137.5 


434.2 


129 


140 



" Unsealed harmonic frequencies. 
St. = stretching; bd. — bending. 
Ref^-i— . Normal mode frequencies observed in vapour phase. 
Ref!^^*^^. Normal mode frequencies were observed in benzene. For comparison, the values of the normal modes are 
calculated at the B3LYP/aug-cc-pVDZ level of theory in solution (benzene) to be 2374.2, 2239.4, 1121.6, 542.10, 500.7, 313.1 

and 121.9 cm~^ for ui.... and respectively. 



TABLE VL Comparison of the IR intensities (in Km/mol) of CICCCN calculated at the HF and DFT levels with corresponding 
experimental values. Here A, B, and C stands for 6-311++g(d,p), aug-cc-PVDZ and aug-cc-PVTZ basis sets respectively. 



li 




RHF 






B3LYP 






B3PW91 




MPa 




A*" 




A 


B 


C 


A 


B 


C 


A 


B 


C 


A 






1 


195.05 


190.18 


185.27 


215.44 


200.78 


194.59 


219.57 


201.19 


195.93 


155.84 


vs 


vs 


2 


38.88 


33.34 


40.98 


27.41 


19.47 


26.30 


31.02 


21.94 


28.92 


18.84 


w 


W 


3 


16.26 


15.10 


14.47 


25.40 


22.60 


21.74 


27.12 


23.74 


22.60 


16.14 


S 


S 


4 


5.49 


4.73 


4.77 


8.15 


7.15 


7.01 


8.23 


7.17 


7.01 


8.11 


m 


m 


5 


9.20 


10.70 


9.62 


5.98 


7.91 


6.93 


5.44 


7.94 


6.84 


8.68 


s 


s 


6 


4.54 


7.37 


10.16 


5.22 


4.26 


7.37 


4.88 


4.03 


7.08 


2.36 


m 


m 


7 


4.87 


5.01 


4.91 


3.90 


3.78 


3.54 


3.86 


3.61 


3.38 


1.26 


vw 


vw 


Sum 


274.29 


266.43 


270.18 


299.0 


265.95 


267.48 


291.50 


269.48 


271.76 


211.23 







" Refp2Si2&. vs = very strong, s — strong, m = medium strong; vw = very weak; w — weak. 
>> Rei,^^. For comparison, the IR intensities are calculated at the B3LYP/aug-cc-pVDZ level of theory in solution (benzene) 
to be 428.7, 21.8, 33.6, 10.2, 12.6, 4.2 and 5.4 cm"^ for Ai... and A7 respectively. 



the BlLYP/TZV(3df,3p) and B3LYP/6-311-|-hg(df,pd) 
levels respectively on the structure as well as on the 
B3PW91/6-311-I— |-g(d,p) optimized geometry are well 
within the scale length of experimental uncertainty. En- 
larging basis sets size could improve the calculation ac- 
curacy and satisfactorily reproduced the experimental 
bending mode frequencies without being scaled with uni- 
form scaling factors, though the stretching frequencies 
are a bit off. The rotational invariants have been sat- 
isfactorily exolained at these levels of theory. The sta- 



theory is not a true minimum rather a saddle point of 
order two and thus scan of the potential energy surface 
should be emphasized in oder to determine the nature 
of the saddle point and verify any possible existence of 
conical intersection between the ground and first excited 
state. 
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